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Q_i| Mass spectra of baryons consisting of two heavy (b or c) and one light quarks are 



calculated in the framework of the relativistic quark model. The light quark-heavy 
diquark structure of the baryon is assumed. Under this assumption the ground and 
excited states of both the diquark and quark-diquark bound system are considered. 



X 

The quark-diquark potential is constructed. The light quark is treated completely 
relativistically, while the expansion in the inverse heavy quark mass is used revealing 
the close similarity with the mass spectra of B and D mesons. We find that the 
relativistic treatment of the light quark plays an important role. The level inversion 
of the p-wave excitations of the light quark in doubly heavy baryons is discussed. 

PACS numbers: 14.20.Lq, 12.39.Ki, 14.20.Mr 



I. INTRODUCTION 



The description of baryons within the constituent quark model is a very important prob- 
lem in quantum chromodynamics (QCD). Since the baryon is a three-body system |], its 
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theory is much more complicated compared to the two-body meson system. Up till now it 
is not even clear which of two main QCD models, Y law or A law, correctly describe the 
nonperturbative (long-range) part of the quark interaction in the baryon [0, |3|. The popular 
quark-diquark picture of a baryon is not universal and does not work in all cases M. The 
success of the heavy quark effective theory (HQET) |J in predicting some properties of the 
heavy-light qQ mesons (B and D) suggests to apply these methods to heavy-light baryons, 
too. The simplest baryonic systems of this kind are the so-called doubly heavy baryons 
(qQQ) i) 0? H H; OH HH) HH • The two heavy quarks (b or c) compose in this case a bound 
diquark system in the antitriplet colour state which serves as a localized colour source. The 
light quark is orbiting around this heavy source at a distance much larger (~ l/w g ) than 
the source size (~ 2 /tuq). Thus the doubly heavy baryons look effectively like a two-body 



bound system and strongly resemble the heavy-light B and D mesons @, 13]. Then the 
HQET expansion in the inverse heavy quark mass can be used. The main distinction of the 
qQQ baryon from the qQ meson is that the QQ colour source though being almost localized 
still is a composite system bearing integer spin values (0, 1, . . .). Hence it follows that the 
interaction of the heavy diquark with the light quark is not point-like but is smeared by the 
form factor expressed through the overlap of the diquark wave functions. Besides this the 
diquark excitations contribute to the baryon excited states. 

In previous approaches for the calculation of doubly heavy baryon masses the expansion 
in inverse powers not only of the heavy quark mass (toq) but also of the light quark mass 
(m q ) is carried out. The estimates of the light quark velocity in these baryons show that the 
light quark is highly relativistic (v/c ~ 0.7 4- 0.8). Thus the nonrelativistic approximation 
is not adequate for the light quark. Here we present a consistent treatment of mass spectra 
of the doubly heavy baryons in the framework of the relativistic quark model based on the 
quasipotential wave equation without employing the expansion in l/m q . Thus the light 
quark is treated fully relativistically. Concerning the heavy diquark (quark) we apply the 
expansion in 1/Mqq (1/wiq). We used a similar approach for the calculation of the mass 



spectra of B and D mesons [14 



The paper is organized as follows. In Sec. II we describe our relativistic quark model 
giving special emphasis to the construction of the quark-quark interaction potential in the 
diquark and the quark-diquark interaction potential in the baryon. In Sec. [ITJ we apply 
our model to the investigation of the heavy diquark properties. The cc and bb diquark 
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mass spectra are calculated. We also determine the diquark interaction vertex with the 
gluon using the quasipotential approach and calculated diquark wave functions. Thus we 
take into account the internal structure of the diquark which considerably modifies the 
quark-diquark potential at small distances and removes fictitious singularities. In Sec. [IV] 
we construct the quasipotential of the interaction of a light quark with a heavy diquark. 
The light quark is treated fully relativistically. We use the expansion in inverse powers of 
the heavy diquark mass to simplify the construction. First we consider the infinitely heavy 
diquark limit and then include the 1/Mqq corrections. In Sec. ^we present our predictions 
for the mass spectra of the ground and excited states of H cc , Ebb, ^ cc and Qbb baryons. We 
consider both the excitations of the light quark and the heavy diquark. The mixing between 
excited baryon states with the same total angular momentum and parity is discussed. For 
H c {, and Q c b baryons, composed from heavy quarks of different flavours we give predictions 
only for ground states, since the excited states of the cb diquark are unstable under the 



emission of soft gluons ||11|| . Moreover, a detailed comparison of our predictions with other 
approaches is given. We reveal the close similarity of the excitations of the light quark in 
a doubly heavy baryon and a heavy-light meson. We also test the fulfillment of different 
relations between mass splittings of doubly heavy baryons with two c or b quarks as well 
as the relations between splittings in the doubly heavy baryons and heavy-light mesons 
following from the heavy quark symmetry. Section |VT| contains our conclusions. 

II. RELATIVISTIC QUARK MODEL 

In the quasipotential approach and quark-diquark picture of doubly heavy baryons the 
interaction of two heavy quarks in a diquark and the light quark interaction with a heavy di- 
quark in a baryon are described by the diquark wave function (^d) of the bound quark-quark 
state and by the baryon wave function (^b) of the bound quark-diquark state respectively, 
which satisfy the quasipotential equation |0| of the Schrodinger type ]16| 

where the relativistic reduced mass is 

E X E 2 M 4 - (mj - m\f 
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and Ei, E 2 are given by 

M 2 — m 2 + m\ M 2 — m\ + m 2 , 

1= 2M ' 2= 2M ' ^ 

here M = E\ + E 2 is the bound state mass (diquark or baryon), rn\ t2 are the masses of heavy 
quarks (Qi and Q 2 ) which form the diquark or of the heavy diquark (d) and light quark (q) 
which form the doubly heavy baryon (B), and p is their relative momentum. In the center 
of mass system the relative momentum squared on mass shell reads 

h 2 (M) _ W 2 - K + m 2 ) 2 }[M 2 - (mi - m 2 ) 2 ] 
1 J 4M 2 ' 1 ' 

The kernel V(p, q; M) in Eq. (JJ) is the quasipotential operator of the quark-quark or 
quark-diquark interaction. It is constructed with the help of the off-mass-shell scattering 
amplitude, projected onto the positive energy states. In the following analysis we closely 
follow the similar construction of the quark-antiquark interaction in mesons which were ex- 
tensively studied in our relativistic quark model |T3|, [IT]]. For the quark-quark interaction 
in a diquark we use the relation Vqq = Vqq/2 arising under the assumption about the 
octet structure of the interaction from the difference of the QQ and QQ colour states. An 
important role in this construction is played by the Lorentz-structure of the confining inter- 
action. In our analysis of mesons while constructing the quasipotential of quark-antiquark 
interaction we adopted that the effective interaction is the sum of the usual one-gluon ex- 
change term with the mixture of long-range vector and scalar linear confining potentials, 
where the vector confining potential contains the Pauli terms. We use the same conventions 
for the construction of the quark-quark and quark-diquark interactions in the baryon. The 
quasipotential is then defined by JT0|, [L7j 

(a) for the quark-quark (QQ) interaction 

V(p, q; M) = Ui(p)u 2 (-p)V(p, q; M) Ul (q)u 2 (-q), (5) 

with 

V(p, q; M) = |a 5 ^(k)^ + (k)r?(k)r^(-k) + ^(k), 

(b) for quark-diquark (qd) interaction 

V(p >q; M) = ^MMu q {p)\a s D,MYu q {q) 

+rAP)uMW(WL { (Vu q (q)MQ) 
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(6) 



where as is the QCD coupling constant, the colour factor is equal to 2/3 for quark- 
quark and 4/3 for quark- diquark interactions, (d(P)\J^\d(Q)) is the vertex the diquark- 
gluon interaction which is discussed in detail below P = (Ed, — p) and Q = (Ed, — q) , 
Ed = (M 2 — m 2 + M|)/(2M) . D^ v is the gluon propagator in the Coulomb gauge 



D 



Qi 



D 



tO 



0. 



(7) 



and k = p — q; 7 M and u(p) are the Dirac matrices and spinors 

u x (p) 



e(p) + m 


( i \ 


\ Hp) 


a- ■ p 


V e(p) +m J 



X 



(8) 



with e(p) = y/p 2 + m' 2 . 

The diquark state in the confining part of the quark-diquark quasipotential @ is de- 
scribed by the wave functions 



^d(p) 




for scalar diquark 

for (axial) vector diquark 



(9) 



where the four vector 



Igrf • PJ 



y£d ■ P)P 



M d (E d (p) + M d 



£d(p) - P = 0, 



(10) 



is the polarization vector of the (axial) vector diquark with momentum p, E d (p) = J p 2 + Mj 
and Ed(0) = (0,£d) is the polarization vector in the diquark rest frame. The effective long- 
range vector vertex of the diquark can be presented in the form 

(P + Q)fj, — scalar diquark 

(11) 



J, 



d',fi 



2^E d (p)E d (q) 

(P + Q)ji _ jj+dyvl 



for scalar diquark 
for (axial) vector diquark 



where k = (0, k) and we neglected the contribution of the chromoquadrupole moment of 
the (axial) vector diquark, which is suppressed by an additional power of k/Md- Here the 
antisymmetric tensor 



-'pa) 



(12) 
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and the (axial) vector diquark spin is given by (Sd-k)u = —isku- We choose the total 
chromomagnetic moment of the (axial) vector diquark fi^ = 2 [HJ |18|| . 
The effective long-range vector vertex of the quark is defined by [14], [17 



T^(k) = 7m + |^v^, fc = (0,k), (13) 

where k is the Pauli interaction constant characterizing the anomalous chromomagnetic 
moment of quarks. In the configuration space the vector and scalar confining potentials in 
the nonrelativistic limit reduce to 

VZnf(r) = (l-OKo„f(r), 

KonfM = eWr), ( 14 ) 

with 

KoBf (r) = KonfM + Ko» = Ar + B, (15) 

where e is the mixing coefficient. 

All the parameters of our model like quark masses, parameters of linear confining potential 
A and B, mixing coefficient e and anomalous chromomagnetic quark moment k were fixed 
from the analysis of heavy quarkonium masses [|l7j and radiative decays |19 |. The quark 



masses m& = 4.88 GeV, m c = 1.55 GeV, m s = 0.50 GeV, m U)d = 0.33 GeV and the 
parameters of the linear potential A = 0.18 GeV 2 and B = —0.30 GeV have standard values 
of quark models. The value of the mixing coefficient of vector and scalar confining potentials 
e = — 1 has been determined from the consideration of the heavy quark expansion |2(| and 
meson radiative decays ||19|| . Finally, the universal Pauli interaction constant k = — 1 has 
been fixed from the analysis of the fine splitting of heavy quarkonia 3 Pj- states |T7| and also 
from the heavy quark expansion [|20|| . Note that the long-range magnetic contribution to the 
potential in our model is proportional to (1 + k) and thus vanishes for the chosen value of 
k = -1. 

III. DIQUARKS IN THE RELATIVISTIC QUARK MODEL 

The quark-quark interaction in the diquark consists of the sum of the spin-independent 
and spin-dependent parts 

V QQ (r) = V^ Q (r) + V^(r). (16) 
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The spin-independent part with the account of v 2 /c 2 corrections including retardation effects 



2f| is given by 



VZUr) 



3 r 
f 



f 
2 
2a 



+ -(Ar + 5) + 





+^W 








tot 



2mi?Ti2 [ 3 r 

2 (P • r 



P 2 + ^ 



+ 



3 r 
1-e 



^1 -e)(l + 2k)At 



£ 1 1 
+ 



x < Ar 



P 



f 

+ 2 



H _ L 2mim2 4 \m 2 m 2 
1-e e ( 1 f N 
2m\m2 4 I mf m 2 



5p 2 



(17) 



where {. . denotes the Weyl ordering of operators and 

Air 



a 3 (fji 



A) = ll - gW/, 



(18) 



/3bln(// 2 /A 2 )' 

with /i fixed to be equal to the reduced mass, rif is a number of flavours and A = 85 MeV. 

The spin- dependent part of the quark-quark potential can be presented in our model |T7 
as follows: 

3 



V$R(r) = a L • S + b 



(S 1 -r)(S 2 -r)-(S 1 -S 2 ) 



+ cS 1 .S 2 + dL.(S 1 -S 5 



(19) 



a = 


vti\mi 1 


b = 


1 


3mim 2 
2 


c = 




3m!m 2 


d = 


1 { 


mim 2 | 



' m\ + m 2 \ 2a s (/i 2 ) lm\+m 2 >A 1 
1 + — ; ^ 3 ^— ; h ~(1 + «) 



4mim 2 / 3 



2a s (fi 2 



+ -(1 + k) 2 (1 -e 



2 4mim 2 r 



(mi + mj) 2 ^ _ A" 



2mim 2 



!^5-(r) + 1(1 + «)"(!- e) f 



m; — m 



2« s (/i 2 



1 A 
27 



m 



2mim 2 



r 



(20) 



where L is the orbital momentum and Si )2 , S = Si + S 2 are the spin momenta. 

Now we can calculate the mass spectra of heavy diquarks with the account of all relativis- 
ts corrections (including retardation effects) of order v 2 /c 2 . For this purpose we substitute 
the quasipotential which is a sum of the spin-independent ( |TTD and spin-dependent flT9"p 
parts into the quasipotential equation (|ip. Then we multiply the resulting expression from 
the left by the quasipotential wave function of a bound state and integrate with respect to 
the relative momentum. Taking into account the accuracy of the calculations, we can use 
for the resulting matrix elements the wave functions of Eq. (||) with the static potential 



l/NR, 
V QQ 1 



r) 



2 a s (/i 2 ) 

3 r 



+ -(Ar + B). 



(21) 
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As a result we obtain the mass formula 

3 



b 2 (M) 

= W + (a)(L-S) + (b)( 



r 



(S 1 -r)(S 2 -r)-(S 1 -S 2 



) + (c)(S 1 -S 2 ) 



+(d) (L-(S 1 -S 2 )), (22) 
where 

w = aa> + g, 

(L • S) = i( J(J + 1) - L(L + 1) - 5(5 + 1)), 

= 6((L ■ S)) 2 + 3(L ■ S) - 2g(S + 1)L(L + 1) 



4(S 1 t)(S 2 t)-(S 1 .S 



2(2L - l)(2L + 3) 
(Si-S 2 ) = i(s(S + l)-^, s = s 1 + s 2 , 

and (a), (6), (c), (d) are the appropriate averages over radial wave functions of Eq. (^ 



is 



We use the notations for the heavy diquark classification: n 2S+l Lj, where n = 1,2, 
a radial quantum number, L is the angular momentum, S = 0, 1 is the total spin of two 
heavy quarks, and J = L — S, L, L + S is the total angular momentum (J = L + S), which 
is considered as the spin of diquark (S^) in the following section. The first term on the 
right-hand side of the mass formula ( p2"D contains all spin-independent contributions, the 
second and the last terms describe the spin-orbit interaction, the third term is responsible 
for the tensor interaction, while the forth term gives the spin-spin interaction. The results 
of our numerical calculations of the mass spectra of cc and bb diquarks are presented in 
Tables | and [IT]. The mass of the ground state of the be diquark in axial vector (l 3 Si) state 
is 



and in scalar (l 1 ^) state is 



M£ = 6.526 GeV 



Ml = 6.519 GeV. 



L bc 

In order to determine the diquark interaction with the gluon field, which takes into ac- 
count the diquark structure, it is necessary to calculate the corresponding matrix element of 
the quark current between diquark states. This diagonal matrix element can be parametrized 
by the following set of elastic form factors 

(a) scalar diquark (S) 

(S(P)\J fl \S(Q)) = h + (k 2 )(P + Q)„ (23) 
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TABLE I: Mass spectrum and mean squared radii of the cc diquark. 



state 


Mass 


/„2\l/2 


btate 


TV /T _ _ _ 

Mass 


/~2\l/2 




(GeV) 


(fm) 




(GeV) 


(fm) 


l 3 Si 


3.226 


0.56 


1^1 


3.460 


0.82 


2 3 ^ 


3.535 


1.02 


2 1 i>i 


3.712 


1.22 


3 3 ^ 


3.782 


1.37 


3 x Pi 


3.928 


1.54 



TABLE II: Mass spectrum and mean squared radii of bb diquark. 



State 


Mass 


(7-2)1/2 


State 


Mass 


( r 2)l/2 




(GeV) 


(fm) 




(GeV) 


(fm) 


l 3 Si 


9.778 


0.37 


l^i 


9.944 


0.57 


2 3 Si 


10.015 


0.71 


2 1 P 1 


10.132 


0.87 


3 3 5i 


10.196 


0.98 


3 1 P 1 


10.305 


1.12 


4 3 5i 


10.369 


1.22 


4 x Pi 


10.453 


1.34 



(b) (axial) vector diquark (V) 

(V(P)\J,\V(Q)) = -[e* d (P)-e d (Q)}hi(k 2 )(P + Q)i 



+h 2 (k 2 ) {[e* d (P) ■ Q}e d ,,(Q) + [e d (Q) ■ P}e* d jP)} 



+h3(k^[e* d (P) ■ Q}[e d {Q) • P](P + Q) M , (24) 

where k = P — Q and £d{P) is the polarization vector of the (axial) vector diquark (|10|). 

In the quasipotential approach, the matrix element of the quark current J M = Q'y^Q 
between the diquark states (d) has the form [E^] 

(d(P)\U0)\d(Q)} = J ^^^(P)r,(p,q)^(q), (25) 

where T At (p,q) is the two-particle vertex function and \l/p are the diquark wave functions 
projected onto the positive energy states of quarks and boosted to the moving reference 
frame with momentum P. The leading contribution to V comes from Fig. |]. Other non- 
leading terms are the consequence of the projection onto the positive-energy states and give 
contributions only at order of 1/rrig for diquarks composed from two heavy quarks and can 
be neglected. Thus we will limit our analysis only to the contribution coming from Fig. [1]. 
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Q2 



Qi 
d 

Q2 



FIG. 1: Lowest order vertex function corresponding to Eq. (|26|). The gluon interaction only 
with one heavy quark is shown. 



The corresponding vertex function is given by 



(26) 



where 22 



rf)(p,q) = n Ql ( Pl )7^Q 1 (gi)(2vr) 3 5(p 2 - q 2 ) + (1 <- 2) 



pi,2 = ei >)^±f>«(py, 

9i,2 = e ll2 (g)^±£n«(Q)g\ 
M d i=1 

and are three four- vectors defined by 

{pi pi pj ~\ 
— , Su + — r , E= VP 2 + M 2 . 
M M(E + M) J 

We substitute the vertex function T^ given by Eq. (|26D in the matrix element (|25|) and 
expand it in I/tuq up to the leading order. Comparing the resulting expressions with the 
form factor decompositions (|23|) and (p4|) we find 



h + {k 2 ) = h 1 (k 2 )=h 2 (k 2 )=2F(k 2 ), 
h(k 2 ) 



0. 



F(k 2 ) 



yEdMd 
E d + M d 



d 3 p 
(2^)3 



*d P + 



2m 



Q-2 



E d + M d 



kU,(p) + (lf,2) 



(27) 



where ^ d = are the di quark wave functions at rest. We calculated corresponding form 
factors F(r)/r which are the Fourier transforms of F(k 2 )/k 2 using the diquark wave func- 
tions found by numerical solving the quasipotential equation. In Fig. |2] the functions F(r) 
for the cc diquark in (IS, IP, 2S, 2P) states are shown as an example. We see that the slope 
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FIG. 2: The form factors F(r) for the cc diquark. The solid curve is for the IS* state, the dashed 
curve for the IP state, the dashed-dotted curve for the 25 state, and the dotted curve for the 2P 
state. 

of F(r) decreases with the increase of the diquark excitation. Our estimates show that this 
form factor can be approximated with a high accuracy by the expression 



F(r) 



-£i — C r 



(28) 



which agrees with previously used approximations |P3| . The values of parameters £ and 



£ for different cc and bb diquark states are given in Tables |I| and [TV]. As we see the 
functions F(r) vanish in the limit r — ► and become unity for large values of r. Such a 
behaviour can be easily understood intuitively. At large distances a diquark can be well 
approximated by a point-like object and its internal structure cannot be resolved. When the 
distance to the diquark decreases the internal structure plays a more important role. As the 
distance approaches zero, the interaction weakens and turns to zero for r = since this point 
coincides with the center of gravity of the two heavy quarks forming the diquark. Thus the 
function F(r) gives an important contribution to the short-range part of the interaction of 
the light quark with the heavy diquark in the baryon and can be neglected for the long-range 
(confining) interaction. It is important to note that the inclusion of such a function removes 
a fictitious singularity 1/r 3 at the origin arising from the one-gluon exchange part of the 
quark-diquark potential when the expansion in inverse powers of the heavy-quark is used. 
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TABLE III: Parameters £ and £ for ground and excited states of cc diquark. 



State 


£ 


c 


State 


£ 


C 




(GeV) 


(GeV 2 ) 




(GeV) 


(GeV 2 ) 


15 


1.30 


0.42 


IP 


0.74 


0.315 


25 


0.67 


0.19 


2P 


0.60 


0.155 


35 


0.57 


0.12 


3P 


0.55 


0.075 




TABLE IV: Parameters £ and £ for 


ground and excited states of bb diquark. 




State 


£ 


C 


State 


£ 


c 




(GeV) 


(GeV 2 ) 




(GeV) 


(GeV 2 ) 


15 


1.30 


1.60 


IP 


0.90 


0.59 


25 


0.85 


0.31 


2P 


0.65 


0.215 


35 


0.66 


0.155 


3P 


0.58 


0.120 


45 


0.56 


0.09 


4P 


0.51 


0.085 



IV. QUASIPOTENTIAL OF THE INTERACTION OF A LIGHT QUARK WITH 

A HEAVY DIQUARK 



The expression for the quasipotential (^j) can, in principle, be used for arbitrary quark and 
diquark masses. The substitution of the Dirac spinors (|8|) and diquark form factors ( P^Q and 
([24]) into (@) results in an extremely nonlocal potential in the configuration space. Clearly, 
it is very hard to deal with such potentials without any simplifying expansion. Fortunately, 
in the case of heavy-diquark light-quark picture of the baryon, one can carry out (following 
HQET) the expansion in inverse powers of the heavy diquark mass The leading terms 
then follow in the limit — > oo. 

A. Infinitely heavy diquark limit 



In the limit — > oo the heavy diquark vertices ( |2"B| ) and ( |2~3| ) have only the zeroth 
component, and the diquark mass and spin decouple from the consideration. As a result we 
get in this limit the quasipotential for the light quark similar to the one in the heavy-light 
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meson in the limit of an infinitely heavy antiquark ||14|| . The only difference consists in the 
extra factor F(k 2 ), defined in fl2"T|), in the one-gluon exchange part which accounts for the 
heavy diquark structure. The quasipotential in this limit is given by 



K(p,q;M) = u q (p)< 



+ Konf(k) 



/V 



7° + 

h 2m, lq 



7,°(7k) 



+ V£ n{ {k)\u q {q). 



(29) 



The resulting interaction is still nonlocal in configuration space. However, taking into 
account that doubly heavy baryons are weakly bound, we can replace e q (p) — > E q = 
(M 2 - Mj + m 2 )/{2M) in the Dirac spinors (§) [|lj] . Such simplifying substitution is widely 



used in quantum electrodynamics [gj, [26[ and introduces only minor corrections of order 
of the ratio of the binding energy (V) to E q . This substitution makes the Fourier transfor- 
mation of the potential fl29|) local. In contrast with the heavy-light meson case no special 
consideration of the one-gluon exchange term is necessary, since the presence of the diquark 
structure described by an extra function F(k 2 ) in Eq. ( |29| ) removes fictitious 1/r 3 singularity 
at the origin in configuration space. 

The resulting local quark-diquark potential for — > oo can be presented in configuration 
space in the following form 

1 



E q + m q 



2E n 



VcouiM + Konf r) + 



(E q + m q ) 2 



p[Vctoui(r) 



conf > 



+ ~ ( ^Coul(r) " KUr) ~ K'onfM 



Eq + mq -A\L f {r)\l 



2m n 



TlV 



conf ' 

L m q 



1 + k)) 

Eg + VTig 



2m n 



1 • s. 



(30) 



where Vcoui(^) = — (4/3)a s F(r)/r is the smeared Coulomb potential. The prime denotes 
differentiation with respect to r, 1 is the orbital momentum, and S q is the spin operator of the 
light quark. Note that the last term in ( |30"D is of the same order as the first two terms and thus 
cannot be treated perturbatively. It is important to note that the quark-diquark potential 
VM d -*oo{r) almost coincides with the quark-antiquark potential in heavy-light (B and D) 



mesons for tuq — > oo [TJ]. The only difference is the presence of the extra factor F{r) in 
Vcoui(^) which accounts for the internal structure of the diquark. This is the consequence of 
the heavy quark (diquark) limit in which its spin and mass decouple from the consideration. 
In the infinitely heavy diquark limit the quasipotential equation ([]]) in configuration space 
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becomes 

(^E^ " *k ) * fl(r) = (r) ^ (r) ' (31) 

and the mass of the baryon is given by M = Md + E q . 

Solving (|3l~l) numerically we get the eigenvalues E q and the baryon wave functions ^b- 
The obtained results are presented in Table [V[ We use the notation ridLn q l(j) for the 
classification of baryon states in the infinitely heavy diquark limit. Here we first give the 
radial quantum number rid and the angular momentum L of the heavy diquark. Then the 
radial quantum number n q , the angular momentum I and the value j of the total angular 
momentum (j = 1 + S q ) of the light quark are shown. We see that the heavy diquark 
spin and mass decouple in the limit Md — > oo, and thus we get the number of degenerated 
states in accord with the heavy quark symmetry prediction. This symmetry predicts also an 
almost equality of corresponding light-quark energies E q for bbq and ccq baryons and their 
nearness in the same limit to the light-quark energies E q of B and D mesons [|l4j . The small 
deviations of the baryon energies from values of the meson energies are connected with the 
different forms of the singularity smearing at r = in the baryon and meson cases. 



B. l/Md corrections 



The heavy quark symmetry degeneracy of states is broken by l/Md corrections. The 
corrections of order l/Md to the potential ( p0[ ) arise from the spatial components of the 
heavy diquark vertex. Other contributions at first order in l/Md come from the one-gluon- 
exchange potential and the vector confining potential, while the scalar potential gives no 
contribution at first order. The resulting l/Md correction to the quark-diquark potential 
(j30|) is given by the following expression 

(a) scalar diquark 



5V, 



1/M d 



E q M d 
1 



P 



A Konf( 



^Coul(r) + Konf( 



P + ^CoulW^: 



V ( r ) | ( 1 + K ) v >v 



conf \ 



IS 



1 I ' 



(32) 



(b) (axial) vector diquark 
6Vl/Md{r) = E^M d 



P 



Coul 



r) + V r 



conf v 



Coull 



2r 
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TART TP A/- Tl-io Trolnoc r~if TP 

lAoijii/ v. ine values 01 H/q 


in the limit — 


-» oo [in \j& v ). 




Baryon 


ccq 


CCS 


bbq 


bbs 


State 


E q 


E s 


E q 


E s 


1S1s(1/2) 


0.491 


0.638 


0.492 


0.641 


lSIp(3/2) 


0.788 


0.906 


0.785 


0.904 


lSIp(l/2) 


0.877 


0.968 


0.880 


0.969 


152s(l/2) 


0.987 


1.080 


0.993 


1.084 


1P1s(1/2) 


0.484 


0.633 


0.489 


0.636 


lPlp(3/2) 


0.793 


0.909 


0.789 


0.906 


lPlp(l/2) 


0.873 


0.965 


0.876 


0.967 


lP2s(l/2) 


0.980 


1.075 


0.984 


1.078 


2Sls(l/2) 


0.481 


0.631 


0.486 


0.634 


2SIp(3/2) 


0.794 


0.909 


0.791 


0.908 


lb Lp\L/ A ) 


U.o 1 1 


U.yoo 


n QV/i 
U.oi4 


u.yoo 


2S2s(l/2) 


0.979 


1.074 


0.982 


1.076 


2Pls(l/2) 


0.479 


0.630 


0.481 


0.631 


351s(l/2) 


0.478 


0.630 


0.480 


0.630 



+- 



i ri 



(1 + k) 



■V, 



IV i 
conf \ 



1S„ 



+ 



+ 3 



l \ r r 

\fyLJr)-VLJT) + (! + *) 
3 

-S g ■ S d + — (S q ■ r)(Srf ■ r) 



is d 
1 



XL(r) - V mn{[ r 



rllV 



']) 



AV Coul (r) + (1 + k)AV c 



conf v 



• S r , 



(33) 



where S = S q + Sd is the total spin, is the diquark spin (which is equal to the total 
angular momentum J of two heavy quarks forming the diquark). The first three terms in 
(|33|) represent spin-independent corrections, the fourth and the fifth terms are responsible 
for the spin-orbit interaction, the sixth one is the tensor interaction and the last one is the 
spin-spin interaction. It is necessary to note that the confining vector interaction gives a 
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contribution to the spin-dependent part which is proportional to (1+k). Thus it vanishes for 
the chosen value of k — — 1, while the confining vector contribution to the spin-independent 
part is nonzero. 

In order to estimate the matrix elements of spin-dependent terms in the 1/Md corrections 
to the quark-diquark potential (|32| ) and ( |33|) as well as different mixings of baryon states, it 
is convenient to use the following relations 

\J;j) = Y,(-l) J+l+Sd+Sq y/&S+l)(2j + l) 

s 

and 

\J;j) = J2(-lY +l+Sd+S ^(2J d + l)(2j + l) 

J d 

where J = j + Sd is the baryon total angular momentum, j = 1 + S q is the light quark total 
angular momentum, S = S q + Sd is the baryon total spin, and 3d = L + S^. 



ltji} |J ' s> (34) 



V. RESULTS AND DISCUSSION 

For the description of the quantum numbers of baryons we use the notations (ridLn q l)J p , 
where we first show the radial quantum number of the diquark {rid = 1, 2, 3 . . .) and its orbital 
momentum by a capital letter (L = S, P, D . . .) , then the radial quantum number of the 
light quark {n q = 1, 2, 3 . . .) and its orbital momentum by a lowercase letter (I — s,p,d . . .), 
and at the end the total angular momentum J and parity P of the baryon. 

The presence of the spin-orbit interaction proportional to 1 • Sd and of the tensor inter- 
action in the quark-diquark potential at 1/Md order (^) results in a mixing of states which 
have the same total angular momentum J and parity but different light quark total momen- 
tum j. For example, the baryon states with diquark in the ground state and light quark 
in the p-wave (I Sip) for J = 1/2 or 3/2 have different values of the light quark angular 
momentum j = 1/2 and 3/2, which mix between themselves. In the case of the ccq baryon 
we have the mixing matrix for J = 1/2 

-55.6 -7.3 \ 

MeV, (36) 

I -8.5 -37.9 J 
with the following eigenvectors 

\(lS2p)l/2'-) = -0.334|j = 3/2) + 0.943|j = 1/2), 
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|(152p)l/2~) = 0.925|j = 3/2) + 0.380|j = 1/2). (37) 
For the ccq baryon with J = 3/2 the mixing matrix is given by 



^ -23.0 18.1 



21.3 18.9 



MeV, (38) 



and the eigenvectors are equal to 



|(l52p)3/2 ") = 0.343|j = 3/2) + 0.939|j = 1/2), 
|(152p)3/2~) = 0.919|j = 3/2) - 0.394|j = 1/2). (39) 

For the bbq baryon we get the mixing matrix for J = 1/2 

MeV, (40) 



-18.0 -2.4 
-2.8 -12.6 



which has eigenvalues 



|(152p)l/2 ") = 0.349|j = 3/2) + 0.937|j = 1/2), 

|(152p)l/2-> = 0.915|j = 3/2) + 0.402|j = 1/2), (41) 



and for J = 3/2 the mixing matrix is 

-7.4 5.q\ 

MeV, (42) 



7.1 6.3 



so that 



|(LS»3/2 ) = 0.341 1 j = 3/2) + 0.940|j = 1/2), 

|(152p)3/2~) = 0.917|j = 3/2) + 0.400|j = 1/2). (43) 

The quasipotential with l/mg corrections is given by the sum of V r mq ^ 00 (r) from (^0[) 
and 5Vi/ mq {r) from (^) and (|33|) . By substituting it in the quasipotential equation ([!]) and 
treating the I/tuq correction term 8V\/ mQ {r) using perturbation theory, we are now able 
to calculate the mass spectra of H cc , E bb , E cb , Q cc , Q bb , Q cb baryons with the account of 



1/rriQ corrections. In Tables [VT|-[rX| we present mass spectra of ground and excited states 
of doubly heavy baryons containing both heavy quarks of the same flavour (c and b) . The 
corresponding level orderings are schematically shown in Figs. ||-|6|. In these figures we first 
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show our predictions for doubly- heavy baryon spectra in the limit when all 1 / corrections 
are neglected [denoting baryon states by ndLn q l(j)]. We see that in this limit the p-wave 
excitations of the light quark are inverted. This means that the mass of the state with 
higher angular momentum j = 3/2 is smaller than the mass of the state with lower angular 
momentum j = 1/2 |13], [14], pTfl . The similar p-level inversion was found previously in 
the mass spectra of heavy-light mesons in the infinitely heavy quark limit [|14|| . Note that 
the pattern of levels of the light quark and level separation in doubly heavy baryons and 
heavy- light mesons almost coincide in these limits. Next we switch on 1/M^ corrections. 
This results in splitting of the degenerate states and mixing of states with different j, which 
have the same total angular momentum J and parity, as it was discussed above. Since the 
diquark has spin one, the states with j = 1/2 split into two different states with J = 1/2 or 
3/2, while the states with j = 3/2 split into three different states with J = 1/2, 3/2 or 5/2. 
The fine splitting between p-levels turns out to be of the same order of magnitude as the gap 
between j = 1/2 and j = 3/2 degenerate multiplets in the infinitely heavy diquark limit. 
The inclusion of 1/M^ corrections leads also to the relative shifts of the baryon levels further 
decreasing this gap. As a result, some of the p-levels from different (initially degenerate) 
multiplets overlap; however, the heavy diquark spin averaged centers remain inverted. The 
resulting picture for the diquark in the ground state is very similar to the one for heavy-light 
mesons []T^| . The purely inverted pattern of p-levels is observed only for the B meson and 
flbb baryons, while in other heavy-light mesons (D, D s , B s ) and doubly heavy baryons 
(5 CC , Q cc ) p-levels from different j multiplets overlap. The absence of the p-level overlap for 
the Qbb baryon in contrast to the B s meson (where we predict a very small overlap of these 
levels [^4|) is explained by the fact that the ratio m s /M^ b is approximately two times smaller 
than m s /rrib and thus it is of order m q /mb. As it was argued in |14[], these ratios determine 
the applicability of the heavy quark limit. 

In Tables [\T| and |VII| we compare our predictions for the ground and excited state masses 
of S cc and Ebb baryons with the predictions of Ref. [fl~l]j . As we see from these Tables our 
predictions are approximately 50 — 150 MeV higher than the estimates of Ref. [PH . One of 
the reasons for the difference between these two predictions for the masses of S cc baryons 
(which is the largest) is the difference in the c quark masses. The mass of the c quark 
in [|TT| is determined form fitting the charmonium spectrum in the quark model where all 
spin-independent relativistic corrections were ignored. However, our estimates show that 
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FIG. 3: Masses of H cc baryons (in GeV). The horizontal dashed line shows the A C D threshold. 
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FIG. 4: Masses of E&& baryons (in GeV). The horizontal dashed line shows the Aj,D threshold. 
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FIG. 5: Masses of O cc baryons (in GeV). The horizontal dashed line shows the K C D S threshold. 
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FIG. 6: Masses of Qbb baryons (in GeV). The horizontal dashed line shows the A),D S threshold. 
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TABLE VI: Mass spectrum of H cc baryons (in GeV). 



State 



(ndLnql)J 1 



Mass 



our 







State 



{n d Ln q l)J } 



Mass 



our 



m 



(lSls)l/2+ 
(LSls)3/2+ 
(151p)l/2- 
(151p)3/2- 
(15lp)l/2'- 
(151p)5/2- 
(151p)3/2'- 



3.620 
3.727 
4.053 
4.101 
4.136 
4.155 
4.196 



3.478 
3.61 
3.927 
4.039 
4.052 
4.047 
4.034 



(LPls)l/2- 
(lPls)3/2~ 
(251s)l/2+ 
(2Sls)3/2+ 
(2Pla)l/2 - 
(2Pl*)3/2 _ 
(351s)l/2+ 



3.838 
3.959 
3.910 
4.027 
4.085 
4.197 
4.154 



3.702 
3.834 
3.812 
3.944 
3.972 
4.104 
4.072 



due to the rather large average value of v 2 /c 2 in charmonium []33| such corrections play an 
important role and give contributions to the charmonium masses of order of 100 MeV. As 
a result the c quark mass found in [|nj is approximately 70 MeV less than in our model. 
For the calculation of the diquark masses we also take into account the spin-independent 
corrections ( |TT| ) to the QQ potential. We find that their contribution is less than in the case 



of charmonium since Vqq = Vqq/2. Thus the cc diquark masses in |TT| are approximately 
50 MeV smaller than in our model. The other main source of the difference is the expansion 
in inverse powers of the light quark mass, which was used in JTTJ] but is not applied in our 
approach, where the light quark is treated fully relativistically. 

In Table |X| we compare our model predictions for the ground state masses of doubly 
heavy baryons with some other predictions || || [IT], |28[ as well as our previous prediction 
||10|| , where the expansion in inverse powers of the heavy and light quark masses was used. 
In general we find a reasonable agreement within 100 MeV between different predictions 
|28|| for the ground state masses of the doubly heavy baryons. The main 



11 



advantage of our present approach is the completely relativistic treatment of the light quark 
and account for the nonlocal composite structure of the diquark. 

For the H c b and Q c t, baryons containing heavy quarks of different flavours (c and b) we 



calculate only the ground state masses. As it was argued in Ref. ||11|| , the excited states 
of heavy diquarks composed of the quarks with different flavours are unstable under the 
emission of soft gluons, and thus the calculation of the excited baryon (cbq and cbs) masses 



22 



TABLE VII: Mass spectrum of H&& baryons (in GeV). 



State 




Mass 


State 


Mass 




(ririLn n l)J p 


our 


U___jJ 


(rLrlLririDJ^ 

\ LL (J I 


our 


H 

ii ii 


(151s)l/2+ 


10.202 


10.093 


(251s)l/2+ 


10.441 


10.373 


(151s)3/2+ 


10.237 


10.133 


(251s)3/2+ 


10.482 


10.413 


(1S1v)1/2- 


10.632 


10.541 


(2Slv)ll2~ 


10.873 




(151»)3/2~ 


10.647 


10.567 


(2Slv)3/2~ 


10.888 




(151p)5/2~ 


10.661 


10.580 


(251p)l/2 /_ 


10.902 




(\S\ri\ll2'~ 


10.675 


10.578 




10.905 




(lSlv)3/2'~ 


10.694 


10.581 


(2Slv)3/2'~ 


10.920 




(152s)l/2+ 


10.832 




(2Pls)l/2~ 


10.563 


10.493 


(152s)3/2+ 


10.860 




(2Pls)3/2 _ 


10.607 


10.533 


(1P1s)1/2~ 


10.368 


10.310 


(351s)l/2+ 


10.630 


10.563 


(lPls)3/2~ 


10.408 


10.343 


(351s)3/2+ 


10.673 




(lPlp)l/2+ 


10.763 




(3Pls)l/2- 


10.744 




(lPlp)3/2+ 


10.779 




(3Pls)3/2- 


10.788 




(lPlp)5/2+ 


10.786 




(451s)l/2+ 


10.812 




(lPlp)l/2'+ 


10.838 




(451s)3/2+ 


10.856 




(lPlp)3/2'+ 


10.856 




(4Pls)l/2~ 


10.900 






TABLE VIII: Mass spectrum 


of Vt cc baryons (in GeV). 




(n r iLn„l)J^' 




Mass 


(n d Ln q l)J p 




Mass 


(151s) 1/2+ 




3.778 


(1P1s)1/2~ 




4.002 


(151s)3/2+ 




3.872 


(lPl«)3/2- 




4.102 


(151p)l/2~ 




4.208 


(251s)l/2+ 




4.075 


(151p)3/2~ 




4.252 


(251s)3/2+ 




4.174 


(15lp)l/2'- 




4.271 


(2Pl«)l/2~ 




4.251 


(l5lp)5/2~ 




4.303 


(2Pls)3/2~ 




4.345 


(lSlp)3/2'- 




4.325 


(351a) 1/2+ 




4.321 
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TABLE IX: Mass spectrum of baryons (in GeV). 



\ ' CI ' Q / 


Mass 


\ ' V U J - J ' V Q V / 


Mass 


(1S1s)1/2+ 


10.359 


(251s)l/2+ 


10.610 


(lSls)3/2+ 


10.389 


(251s)3/2+ 


10.645 


(lSlv)l/2~ 


10.771 


(251p)l/2 _ 


11.011 


(lSlp)3/2 _ 


10.785 


(2Slp)3/2~ 


11.025 


(lSlp)5/2~ 


10.798 


(2Slp)l/2 /_ 


11.035 


dSl»H/2 /_ 


10.804 


(2Sl?/)5/2 _ 


11.040 


(lSlp)3/2 /_ 


10.821 


(2Slp)3/2 /_ 


11.051 


(15'2s)l/2 + 


10.970 


(2Pls)l/2~ 


10.738 


(lS2s)3/2+ 


10.992 


(2Pls)3/2~ 


10.775 


(lPls)l/2 _ 


10.532 


(3Sls)l/2+ 


10.806 


(lPls)3/2 _ 


10.566 


(3Sls)3/2+ 


10.843 


(lPlp)l/2+ 


10.914 


(3Pla)l/2 _ 


10.924 


(lPlp)3/2+ 


10.928 


(3Pls)3/2- 


10.961 


(lPlp)5/2+ 


10.937 


(4Sls)l/2+ 


10.994 


(lPlp)l/2'+ 


10.971 


(4Sls)3/2+ 


11.031 


(lPlp)3/2'+ 


10.986 


(4Pla)l/2~ 


11.083 



is not justified in the quark- diquark scheme. We get the following predictions for the masses 
of the ground state cbq baryons: 

(l<Sls)l/2 + states with the axial vector and scalar cb diquarks respectively 

M(E cb ) = 6.933 GeV, M(E' cb ) = 6.963 GeV, 

(LSls)3/2+ state 

M(E* cb ) = 6.980 GeV, 

and for cbs baryons: 

(lSls)l/2 + states with the axial vector and scalar cb diquarks 

M(n cb ) = 7.088 GeV, M(Q' cb ) = 7.116 GeV, 
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TABLE X: Mass spectrum of ground states of doubly heavy baryons (in GeV). Comparison of 
different predictions. {QQ} denotes the diquark in the axial vector state and [QQ] denotes diquark 
in the scalar state. 



Baryon 


Quark 
content 


J p 


Present 
work 





E3 


| 

1— 1 


i 


m 


1 — 'CC 


{cc}q 


1/2+ 


3.620 


3.478 


3.66 


3.66 


3.61 


3.69 


1 — CC 


{cc}q 


3/2+ 


3.727 


3.61 


3.81 


3.74 


3.68 




^cc 


{cc}s 


1/2+ 


3.778 


3.59 


3.76 


3.74 


3.71 


3.86 


Q* 

CC 


{cc}s 


3/2+ 


3.872 


3.69 


3.89 


3.82 


3.76 




^bb 


{bb}q 


1/2+ 


10.202 


10.093 


10.23 


10.34 




10.16 


C* 
^bb 


{bb}q 


3/2+ 


10.237 


10.133 


10.28 


10.37 






&>bb 


{bb}s 


1/2+ 


10.359 


10.18 


10.32 


10.37 




10.34 


O* 
lL bb 


{bb}s 


3/2+ 


10.389 


10.20 


10.36 


10.40 






^cb 


{cb}q 


1/2+ 


6.933 


6.82 


6.95 


7.04 




6.96 


^cb 


[cb]q 


1/2+ 


6.963 


6.85 


7.00 


6.99 






^cb 


{cb}q 


3/2+ 


6.980 


6.90 


7.02 


7.06 






^cb 


{cb}s 


1/2+ 


7.088 


6.91 


7.05 


7.09 




7.13 


Kb 


[cb]s 


1/2+ 


7.116 


6.93 


7.09 


7.06 






Kb 


{cb}s 


3/2+ 


7.130 


6.99 


7.11 


7.12 







(lSls)3/2+ state 

M(Vl* cb ) = 7.130 GeV. 

Now we compare our results with the model-independent predictions of the heavy quark 
effective theory. The heavy quark symmetry predicts simple relations between the spin 
averaged masses of doubly heavy baryons with the accuracy of order 1/M^ 

AM 1>2 = M 1 (E bb ) - Mi (E,*) = M 2 (E bb ) - M 2 (E CC ) = M l {yi hh ) - M 1 (Q CC ) 

= M 2 (n bb ) - M 2 (Sl cc ) = Mi - M^ c = AM d , (44) 

where the spin-averaged masses are 

Mi = (Mi/a + 2M 3/2 )/3, 



25 



TABLE XI: Differences between spin-averaged masses of doubly heavy baryons denned in Eq. (pPfl ) 
(in GeV). 



Baryon state 


AMi(S) 


AM 2 (H) 


AMi(O) 


AM 2 (fi) 


AM d 


1SU 


6.534 




6.538 




6.552 


ISlp 


6.512 


6.532 




6.519 


6.552 


lPls 


6.476 




6.486 




6.484 


251s 


6.480 




6.492 




6.480 



M 2 = (M 1/2 + 2M 3/2 + 3M 5/2 )/6 

and Mj are the masses of baryons with total angular momentum J. Mqq are the masses of 
diquarks in definite states. The numerical results are presented in Table |X1] (only the states 
below threshold are considered). We see that the equalities in Eq. ([P]) are satisfied with 
good accuracy for the baryons with the heavy diquark and light quark both in ground and 
excited states. 

It follows from the heavy quark symmetry that the hyperfme mass splittings of initially 
degenerate light quark states 

AM(Hqq) = M 3/2 (E QQ )-M 1/2 (E QQ ), 

AM(Q QQ ) ee M 3/2 (Q QQ ) - M 1/2 (Qqq), (45) 
should scale with the diquark masses: 

AM(S CC ) = R AM(E bb ), 

AM(n cc ) = R AM(n bb ), (46) 

where R = M^ h /Mf c is the ratio of diquark masses. Our model predictions for these splittings 
are displayed in Table (XTI1 Again we see that heavy quark symmetry relations are satisfied 



with high accuracy. 

The close similarity of the interaction of the light quark with the heavy quark in the 
heavy-light mesons and with the heavy diquark in the doubly heavy baryons produces very 
simple relations between the meson and baryon mass splittings [p], 29, j3(|. In fact for the 



ground state hyperfme splittings of mesons and baryons we obtain adopting the approximate 
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TABLE XII: Hyperfine splittings (in MeV) of the doubly heavy baryons for the states with the 
light quark angular momentum j = 1/2. 



Baryon state 


R 


AM(E bb ) 


R AM(E bb ) 


AM(E CC ) 


AM(Q bb ) 


R AM(Q bb ) 


AM(fl cc ) 


1S1s[3/2 - 1/2] 


3.03 


35 


106 


107 


30 


91 


94 


151p[3/2 - 1/2] 


3.03 


19 


58 


60 


17 


52 


54 


lPls[3/2 - 1/2] 


2.87 


40 


115 


121 


34 


98 


100 


2Sls[3/2 - 1/2] 


2.83 


41 


116 


117 


35 


99 


99 



TABLE XIII: Comparison of hyperfine splittings (in MeV) in doubly heavy baryons and heavy- light 
mesons. Experimental values for the hyperfine splittings in mesons are taken from Ref. pi]]. 



AM(E CC ) 


|AM™ P AM(E bb ) 


|AM^ xp AM(n cc ) 


|AM** P AM(H 66 ) 




107 


106 35 


34 94 


108 30 


35 



relation Mqq = 2m Q 

AM(H QQ ) = l^r^M B , D - \/AM B ^ 

AM(n QQ ) = |aM Bs , Ds) (47) 

where the factor 3/2 is just the ratio of the baryon and meson spin matrix elements. The 
numerical fulfillment of relations (H7|) is shown in Table |Xlli. 



VI. CONCLUSIONS 



In this paper we calculated the masses of the ground and excited states of the doubly 
heavy baryons on the basis of the quark-diquark approximation in the framework of the 
relativistic quark model. The orbital and radial excitations both of the heavy diquark and 
the light quark were considered. The main advantage of the proposed approach consists in 
the fully relativistic treatment of the light quark (u, d, s) dynamics and in the account for the 
internal structure of the diquark in the short-range quark-diquark interaction. We apply only 
the expansion in inverse powers of the heavy diquark mass (M^ 6 ,M^ C ), which considerably 
simplifies calculations. The infinitely heavy diquark limit as well as the first order 1/Md spin- 
independent and spin-dependent contributions were considered. A close similarity between 
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excitations of the light quark in the doubly heavy baryons and heavy-light mesons was 
demonstrated. In the infinitely heavy (di) quark limit the only difference originates from 
the internal structure of the diquark which is important at small distances. The first order 
contributions to the heavy (di) quark expansion explicitly depend on the values of the heavy 
diquark (boson) and heavy quark (fermion) spins and masses (Mqq « 2m Q ). This results 
in the different number of levels to which the initially degenerate states split as well as their 
ordering. Our model respects the constraints imposed by heavy quark symmetry on the 
number of levels and their splittings. 

We find that the p-wave levels of the light quark which correspond to heavy diquark spin 
multiplets with j = 1/2 and j = 3/2 are inverted in the infinitely heavy diquark limit. The 
origin of this inversion is the following. The confining potential contribution to the spin- 
orbit term in (^) exceeds the one-gluon exchange contribution. Thus the sign before the 
spin-orbit term is negative, and the level inversion emerges. However, the corrections, 
which produce the hyperfine splittings of these multiplets, are substantial. As a result the 
purely inverted pattern of p levels for the heavy diquark in the ground state occurs only 
for the doubly heavy baryons "Ebb and Qbb- For S cc and Q cc baryons the levels from these 
multiplets overlap. The similar pattern was previously found in our model for the heavy-light 
mesons ||14|| . 

We plan to use the found wave functions of doubly heavy baryons for the calculation of 
semileptonic and nonleptonic decays to H c b and S c b to S cc baryons. The corresponding 
baryonic Isgur-Wise functions || ^2| will be determined. 
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